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INTRODUCTION 


1 . 1  Summary 

This  report  details  results  of  research  performed  during  the  period  June 
1,  1988  through  December  31,  1988  under  AFOSR  contract  no.  F49620-86-K-0016. 
The  project  was  initiated  in  June  1986  and  covers  a  duration  of  forty  months. 

The  general  objective  of  this  research  is  to  improve  on  existing 
theoretical  models  for  predicting  the  response  of  inelastic  aerospace 
structural  components  subjected  to  hostile  thermal  environments  with  emphasis 
on  transient  temperature  conditions,  radiation  boundary  conditions,  extremely 
rapid  heating  rates,  and  possible  phase  change  of  the  materials  involved. 

1.2  Statement  of  Work 

Experimental  and  theoretical  research  are  being  performed  to  characterize 
the  response  of  structural  components  subjected  to  transient  temperature 
conditions  resulting  in  inelastic  material  behavior.  The  research  is  being 
performed  in  the  following  stages: 

1)  theoretical  development  of  thermodynamic  constraints  on  inelastic  materials 
under  transient  temperature  conditions; 

2)  development  of  modified  heat  conduction  equations  to  account  for  two-way 
thermomechanical  coupling  in  these  inelastic  materials; 

3)  experimentation  to  determine  further  constraints  on  inelastic  materials 
under  transient  temperature  conditions; 


1 


4)  development  of  multi -dimensional  theoretical  algorithms  for  predicting 
response  of  the  inelastic  structural  components  described  above;  and 

5)  experimentation  to  verify  the  theoretical  algorithms  described  in  item  4). 

Items  1)  through  4)  above  are  being  performed  entirely  on  the  main  campus  at 
Texas  A&M  University.  Item  5)  will  be  performed  both  at  the  Air  Force  Wright 
Aeronautics  Laboratory  at  Wright  Patterson  Air  Force  Base  and  Texas  A&M 
Universitj.  Details  of  this  interaction  will  be  described  further  below. 

RESEARCH  COMPLETED 

2.1  Summary  of  Completed  Research 

Research  during  the  period  June  1  through  December  31,  1988  is  summarized 
below: 

1)  additional  computational  results  have  been  obtained  utilizing  the  one¬ 
way  coupled  model  discussed  in  reference  1; 

2)  two  new  two-way  coupled  models  are  under  development  at  this  time; 

3)  the  initial  plate  experiments  have  been  performed  at  AFWAL  and  the 
second  set  of  plate  experiments  is  being  prepared  at  this  time; 


4)  the  constitutive  tests  on  Hastelloy  X  are  underway;  and 


5)  computational  facilities  have  been  upgraded  to  perform  the  theoretical 


analyses. 

2.2  Expected  Research 

The  following  results  are  expected  between  now  and  Sept.  30,  1989: 

1)  the  plate  experiments  will  be  completed  at  AFWAL; 

2)  the  one-way  coupled  model  will  be  compared  to  the  plate  experiments; 

3)  the  two-way  coupled  model  will  be  completed  and  compared  to  the  plate 
experiments;  and 

4)  the  constitutive  tests  on  Hastelloy  will  be  completed. 


2.3  Theoretical  Developments 

A  brief  description  of  theoretical  progress  over  the  past  six  months  is 
given  below. 

2.3.1  One-Way  Coupled  Model 

The  one-way  coupled  model  has  been  completed  and  documented  in  reference 
1  (See  Appendix  6.1).  Recently  this  model  has  been  used  to  predict  the 
response  of  a  circular  plate  subjected  to  the  instantaneous  heat  pulse  shown 
in  Fig.  1.  As  shown  in  Fig.  2,  the  model  predicts  a  dynamic  response  which  is 
not  unlike  resonance.  Figure  3  shows  the  displacement  profile  during  a  single 
cycle  of  response,  demonstrating  that  a  permanent  out-of-plane 
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deformation  caused  by  the  thermal  field  constrains  this  dynamic  response. 
Figurjs  4  through  7  show  that  significant  stresses  occur  which  causes 
substantial  inelastic  strain  near  the  center  of  the  plate. 


In  our  next  series  of  experiments  we  will  be  attempting  to  determine  the 
veracity  of  this  dynamic  response  as  predicted  by  the  model. 

2.3.2  Two-way  Coupled  Models 

We  are  currently  placing  the  emphasis  in  model  development  on  the 
construction  of  two-way  coupled  algorithms  for  analyzing  the  plate  problem. 

The  two-way  coupled  problem  is  constructed  by  first  postulating 
constitutive  equations  to  be  equations  of  state  in  the  strain, 
temperature,  T,  temperature  gradient,  gj^,  and  internal  state,  u[^^,  as  follows 
121: 

the  stress  tensor: 

°ij  =  (1) 

the  internal  energy: 

the  entropy: 

^  ^ ^ ^ k 1 ’ ^ ’ ^k ’^k  1  ^ 
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Fig.  7 


Inelastic  Strain  Profile  for  Circular  Plate  with  Instantaneous  Heat  Input 
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the  heat  flux  vector: 


(4) 


and 


a^j  ^ij^^kl’^’^k’*^kl’*^kl^  ’  ^  l«***»o,  p  l,...,n  (5) 

where  are  a  set  of  n  internal  states  variables.  Now  define  the  Helmholtz 
free  energy 


h  u-Ts  =  h(L|,,  ,  r,g^^,aj^) 

Therefore, 


(6) 


u  3h 

'' "  -jr 


kl 


■kl 


ah  j  3h 

3T  3g. 


3h 


dCi 


*kl 


kl 


(7) 


Substituting  the  above  into  the  first  law  of  thermodynamics  and  this  result 
into  the  second  law  of  thermodynamics  gives 


\ 


12 


oTs, 


Ti-i  V^\  3h  /  -|.  U 

’  '  ‘^ae- j^^kr^*9k’“kl 


-  I  p 


j  h 
3  Q  . 


(k-k]  .T,g 


3h 

aajj 


^‘-'kr^’Sk’^kl^l  «5j(eki 


(qi(Pkr'^.g 


k’^kl 


"  )/T)g^  i  0 


(8) 


Using  the  Coleman-Mizel  |31  procedure,  it  now  follows  that 


lii  =  0  => 


391 


h=h(e^^^,T,ai;T) 


kV 


(9) 


By  a  similar  argument 


s 


-ah 

aT 


S  s(e|^^  ,T,aJ^) 


(10) 


Also, 


1 J 


3h 

’Se. 


IJ 


IJ 


'i 


(11) 


Note  that  although  the  above  is  similar  to  the  result  obtained  for  elastic 
materials,  h  is  no  longer  a  potential  for  the  stress  tensor  because  it  is  path 
dependent  due  to  its  dependence  on  the  internal  state. 

Inequality  (8)  now  reduces  to 


pT, 


ah 


q.g./T  >  0 


(12) 


where  the  first  term  represents  internal  dissipation  and  the  last  term  is  heat 
conduction  dissipation.  Note  that  the  internal  dissipation  cannot  be  set  to 
zero  because  although  the  rate  of  change  of  the  internal  state  variables  may 
be  specified,  the  actual  internal  state  at  any  time  cannot  be  specified. 

It  can  also  be  shown  that  [3] 

di  =  gj  +  H.O.T.  (13) 

Thus,  if  internal  state  variable  growth  laws  (5)  can  be  determined  the 
problem  will  be  completely  specified  by  construction  of  the  Helmholtz  free 
energy  function  (equation  (9)).  v 

Now  consider  a  special  case  of  equation  (9)  which  is  found  to  be  suitable 
for  many  materials.  Let 

h  =  h(c^^,T,.J^)  (14) 
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where  e^i,  called  the  inelastic  strain  tensor,  is  equivalent  to 

Expanding  h  in  a  second  order  Taylor  series  in  its  arguments  gives 


h  =  — |A-*BT'  +  C.  -t  -  .+iiO  .  ..  ,c  -  •£.  1+E.  -e?  .+F.  .. 

p'  ij  1J  *  ijkl  1J  kl  IJ  1J  ijkl  ij  kl 


(15) 


Substituting  the  above  into  (11)  results  in 


0..  =  C**+D.-,  -i+H.-T* 

ij  ij  ijkl  kl  ijkl  kl  ij 


(16) 


The  above  may  be  written  equivalently 


’i  j  "  °i  j^°i jkJ  'kT'^kl  “kl^'''"'''c 


(17) 


where  o  .  is  the  residual  stress,  D...,is  the  linear  elastic  modulus  tensor, 

n  j  I  jK I 

and  is  the  thermal  expansion  coefficient  tensor.  From  (17)  it  is  clear 


that 


c?  •  =  E-  .-e^  .-a-  .(T-T  ) 
IJ  IJ  TJ  TJ^  0^ 


(18) 


where  is  the  elastic  strain  tensor,  given  by 


^  n-l  t  R  i 

^ij  "  ^ijkr°kr°kr 


(19) 
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The  above  equations  will  be  shown  to  be  suitable  for  characterizing  the 
material  behavior  of  elastic-plastic  media  with  internal  state  variables. 

To  obtain  the  coupled  heat  conduction  equation  recall  the  first  law: 


pu  =  G .  .  c  .  .  -  q  .  .  +  pr 
ij  ij 


(20) 


Therefore,  due  to  definition  (6)  we  obtain 


p(ntTs.Ts)  »  +  pf- 


(21) 


From  (9)  we  know  that 


bh  •  b h  Y  .  bh  •  n 


u  O'*  OHt.wII  ’ 

h  •*'  G  .  4  +  T  +  “ '  Ol  .  . 

’‘iJ  a..;,  ’J 

Substituting  (22)  into  (21)  gives 


(22) 


■  ^ij  °  (|t  ^  ^ 


.  ah  •  bh  -n 

"t"  0  ~  ~~  Q  4  D  — —  CX  •  •  Q  • 

a  n  IJ  ^J.J 

1  au.j. 


pr  =  0 


(23) 


Utilizing  (9)  through  (11)  results  in 


T*  ah'P 

„T5  ,  „  -  .  a-  ,  -  „r  .  0. 

’“ij 


(24) 


Utilizing  (3)  and  (10)  gives 
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Note  that  for  an  elastic  material  d"? .  =  0  and  the  above  reduces  to  the 

^  J 

coupled  heat  conduction  equation  for  elastic  media.  Since  terms  1  and  2  could 
be  large  even  under  non-inertiul  conditions,  they  should  be  considered 
carefully  in  inelastic  problems. 

Finally,  substituting  (17)  into  the  above  gives 
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.1 


■  ^ijkl  ■  ‘‘kl  ‘"kl  ”  ^ijkl  “ij  \l  ^ijkl  “ij  ^  ^ 


kl 


+  p  C^  T  -  (k..  g.)^.  -  pr  =  0 


(29) 


The  above  equation  is  the  two-way  coupled  heat  conduction  equation.  Note 
that  a  change  in  either  the  state  of  strain,  or  the  internal  state, 
ajj,  has  an  effect  on  the  temperature,  T. 

The  authors  nave  previously  studied  the  above  equation  for  a  point  mass 
14)  and  in  one  dimension  [51.  The  inclusion  of  two-way  coupling  adds 
considerable  complexity  to  the  computational  algorithm  for  solving  the  problem 
in  any  number  of  dimensions. 

We  are  currently  developing  two  algorithms  to  account  for  this 
thermomechanical  coupling.  In  the  first,  we  have  reduced  all  of  the  field 
equations,  including  (29),  to  two-dimensional  axisymmetric  form  and  are 
utilizing  continuum  elements  to  obtain  a  solution  for  both  the  temperature  and 
displacement  fields. 

In  the  second  method,  we  are  modifying  the  one-way  coupled  code 
previously  developed  by  the  authors  [5].  We  are  using  an  operator  splitting 
method  previously  proposed  by  Oden  [6].  In  this  approach,  the  problem  is 
solved  on  each  time  increment  by  assuming  one-way  coupling  and  the  coupling 
terms  are  then  included  in  the  thermal  analysis  until  iterative  convergence  is 
obtained  on  each  time  step.  We  hope  to  have  one  or  both  of  these  algorithms 
operational  by  June,  1989.  ' 


2.3.3  Computational  Facilities  Development 

It  was  apparent  from  the  inordinate  computational  times  required  to 
operate  the  one-way  coupled  code  that  it  would  be  necessary  to  improve  our 
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computational  efficiency  in  order  to  develop  two-way  coupled  algorithms  15]. 
Therefore,  the  authors  have  spent  the  last  several  months  upgrading  our 
computer  capabilities.  We  have  purchased  a  microVAX  to  be  utilized  as  a 
dedicated  machine  for  use  with  the  algorithms  developed  as  part  of  the  current 
contract.  This  system  has  now  been  connected  to  the  campus  VAX  network,  and 
initial  indications  are  that  the  computational  times  will  be  decreased 
somewhat.  However,  because  the  improved  efficiency  has  not  been  as  much  as  we 
had  hoped  for,  we  have  undertaken  to  streamline  the  one-way  coupled  code 
somewhat  before  proceeding  to  implement  the  two-way  coupling.  These 
efficiency  measures  are  nearing  completion  at  this  time. 


2.4  Experimental  Programs 

The  experimental  program  has  been  divided  into  two  major  categories:  1) 
constitutive  testing  of  uniaxial  specimens  under  steady  state  and  transient 
temperature  conditions;  and  2)  structural  testing  of  plates  subjected  to  rapid 
external  heating.  The  experimental  efforts  are  discussed  in  greater  detail  in 
the  following  two  subsections. 

The  constitutive  tests  are  used  to  evaluate  the  material  parameters  for 
the  thermoviscoplastic  constitutive  models  used  in  the  structural  analysis 
program  described  in  Section  2.2.  Currently,  these  models  will  utilize 
constants  obtained  from  tests  conducted  at  several  elevated  steady  state 
temperatures.  Material  response  at  intermediate  temperatures  is  then  computed 
using  interpolated  values  of  the  parameters.  This  approach  may  lead  to 
erroneous  predictions  if  the  material  undergoes  a  significant  phase  change 
during  a  temperature  transient.  To  investigate  this  phenomenon,  a  set  of 
experiments  is  being  developed  and  performed  which  will  enable  us  to  extend 
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tne  current  constitutive  theories  to  predict  material  behavior  under  transient 
temperature  conditions. 

The  structural  testing  of  plates  subjected  to  rapid  external  heating  is 
being  performed  in  order  to  critically  compare  the  theory  developed  herein. 
Primary  emphasis  is  being  placed  on  prediction  of  the  plate  response  for  both 
temperature  and  displacement  fields.  In  addition,  more  physical  insight  to 
the  problem  we  are  modelling  is  being  obtained. 

2.4.1  Constitutive  Testing 

Uniaxial  constitutive  experiments  are  currently  being  conducted  in  a 
steady  state  temperature  environment  in  order  to  verify  and  supplement  the 
material  parameter  data  base  obtained  from  |7].  Hastelloy  X  has  been  selected 
as  the  candidate  material  for  this  research  primarily  because  of  the 
availability  of  such  information.  Therefore,  upon  completion  of  these  tests, 
material  constants  for  Hastelloy  X  at  room  temperature,  1000°F,  1200'F, 
1600  F,  and  1800  F  will  be  used  to  perform  the  structural  analysis  discussed 
in  the  previous  sections. 

The  uniaxial  specimens  used  for  testing  are  standard  LCF  uniform-gage 
test  section  specimens,  machined  to  ASTM  specifications.  Typical  tests  being 
performed  include:  1)  creep  tests;  2)  stress  relaxation  tests;  3)  monotonic 
constant  strain  rate  tests;  4)  fully  reversed  cyclic  tests;  and  5)  stress  drop 
tests  during  monotonic  and  cyclic  loading.  All  experiments  are  being 
performed  at  Texas  A&M  University  using  an  MTS  880  electro-hydraulic  testing 
machine  and  clam  shell  oven. 

The  three  zone  clainsHell  oven  and  closed  loop  control  system  allows  for 
accurate  temperature  control  in  the  specimen  gage  se:tion  with  nominal  spatial 
temperature  gradients.  Specimen  temperature  is  obtained  using  three,  contact 
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mounted,  K-type  thermocouples.  While  intrinsically  mounted  thermocouples 
would  allow  for  more  accurate  temperature  measurement,  past  experience  has 
shown  that  this  method  of  attachment  severely  degrades  the  fatigue  life  of  the 
spec imen . 

Constitutive  tests  of  uniaxial  specimens  under  transient  temperature 
conditions  is  also  underway,  in  order  to  extend  the  validity  of  current  models 
to  high  thermal  rates.  Specimens  specifically  designed  for  combined 
thennomechan ical  loading  are  being  evaluated  for  this  phase  of  testing.  The 
specimen  geometry  is  that  of  an  LCF  efficiency  button  head  specimen,  machined 
to  ASFM  specifications,  lo  study  the  uniformity  of  the  transverse  temperature 
gradient,  specimens  have  been  fabricated  to  have  either  a  solid  or  hollow 
core.  In  order  to  provide  test  to  test  compatibility,  hastelloy-X  is  also 
being  used  in  this  phase  of  testing. 

Another  difficult  aspect  of  transient  temperature  testing  is  minimization 
of  the  longitudinal  thermal  gradient  during  a  thermal  transient.  A  quad- 
elliptical  quartz  lamp  furnace  is  being  used  and  will  hopefully  alleviate  this 
potential  problem.  Approximately  thirty  K-type  thermocouples  are  being  used 
to  monitor  the  temperature  fields  in  the  gage  section  of  the  specimen. 

In  order  to  obtain  useful  information  about  the  thermomechanical 
constitutive  behavior  of  inelastic  materials  under  transient  temperature 
conditions,  it  is  necessary  to  impose  extremely  complex  strain  histories,  in 
addition  to  the  thermal  histories.  The  MTS  880  electo-hydrau 1 i c  testing 
machine  is  being  modified  to  simultaneously  measure  mechanical  and  thermal 
data  and  control  (in  a  closed  loop  fashion)  both  aspects  of  the  experiment. 
Therefore,  one  of  the  main  goals  of  the  phase  of  testing  is  to  identify 
acceptable  experimental  techniques  for  these  type  of  experiments. 
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2.4.2  Plate  Testing 


Indications  from  the  first  set  of  plate  experiments  performed  at  AFWAL 
are  that  the  acquired  data  are  not  representative  of  the  actual  experimental 
output.  This  is  believed  to  be  due  to  an  inordinate  amount  of  electronic 
noise  which  masked  the  output  produced  by  both  the  LVDT's  and  the 
tnerniocoup  les.  We  hypothesize  at  this  time  that  this  was  caused  by  the 
installation  of  a  second  laser  on  the  day  that  the  experiments  were 
performed.  Therefore,  we  are  currently  planning  to  rerun  the  tests  in  the 
near  future.  At  the  present  time,  the  technicians  onsite  are  preparing  the 
experiments  arid  debugging  the  electronic  data  acquisition  units.  The  second 
set  of  tests  are  scheduled  to  be  performed  in  March. 

2. 5  Conclusions 

There  exists  a  need  to  develop  models  capable  of  predicting  the  response 
of  aerospace  structures  to  complex  thermomechanical  inputs  which  produce 
significant  material  inelasticity  and  resulting  two-way  thermomechanical 
coupling  and  to  verify  these  models  against  controlled  experiments. 

This  report  documents  progress  made  during  the  past  six  months  towards 
achieving  this  goal. 
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ABSTRACT 


A  computational  analysis  is  presented  herein  for  an  axisymmetrir  plate 
subjected  to  rapid  external  heating.  The  analysis  is  complicated  by  several 
forms  of  nonlinearity,  including  radiation  boundary  conditions,  material 
viscoplasticity  and  geometric  nonlinearity.  The  mechanical  constitution  is 
extremely  complc/.,  involving  history  and  rate  dependence  which  causes  the 
constitutive  equations  to  be  mathematically  stiff.  Results  are  obtained  for 
two  viscoplasticity  models  available  in  the  current  literature. 

The  solution  utili2ed  herein  adopts  the  finite  element  method  in  spatial 
coordinates  and  standard  finite  differencing  in  time.  Iterative  techniques 
are  used  to  account  for  nonlinearity.  Results  are  obtained  for  a  circular 
plate  subjected  to  a  high  energy  instantaneous  heat  source  applied 
axisymmetrical ly  to  one  side  of  the  plate.  Sensitivity  studies  are  conducted 
to  determine  the  effects  of  various  heat  intensities  and  plate  thicknesses  on 
temperature,  displacements,  and  stresses.  It  is  found  that  material 
viscoplasticity  and  geometric  nonlinearity  contribute  significantly  to  the 
predicted  response  of  the  plate. 

INTRODUCTION 

Plates  are  often  subjected  to  rapid  heating  capable  of  producing  highly 
nonlinear  structural  response.  Examples  are  structural  components  subjected 
to  laser  heating,  the  skins  of  aerodynamic  vehicles  in  hypersonic  flight,  and 
hot  gas  turbine  engine  components.  Du^  to  the  elevated  temperature 
environment  encountered  in  these  applications  it  is  often  necessary  to  include 
thermal  effects  in  the  structural  analysis. 

several  ti.  oretical  solutions  for  a  heated  thin  plate  have  been  reported 
in  the  literature.  In  most  cases  e-'ther  classical  plate  theory  [1,2]  or  large 
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deflection  plate  theory  [3-6J  was  utilized  to  obtain  linear  thermoelastic 
solutions.  In  severe  thermal  environments,  however,  the  plate  response  may 
become  highly  nonlinear.  Large  variations  in  the  temperature  field 
necessitate  inclusion  of  temperature  dependent  thermal  conductivity. 
Furthermore, the  nonuniformly  distributed  thermal  strain  induces  geometric 
nonlinearity  and  localized  inelastic  stresses.  It  is  thus  a  formidable  task 
to  extend  the  previous  solutions  to  account  for  these  nonl inearities._ 

Recently,  attempts  have  been  made  to  numerically  approximate  the  solution 
for  a  heated  elastic  plate  17,8|  by  the  boundary  element  method.  Kawakami  and 
Shiojiro  19|  attempted  to  implement  AOINAT/ADINA  with  the  capability  of 
performing  heat  transfer  analysis  and  thermoelastic-plastic  analysis  for 
plate/shell  elements.  In  their  research,  plasticity  in  the  element  is  assumed 
only  when  the  cross  section  is  totally  plastic.  Therefore,  their  solution 
does  not  apply  to  a  plate  subjected  to  rapid  high  energy  heating  on  one 
surface. 

The  authori,  have  reported  a  computational  analysis  for  an  axisymmetric 
viscoplastic  plate  subjected  to  rapid  external  heating  [lO].  To  our  knowledge 
this  is  the  first  plate  analysis  involving  both  heat  transfer  and 
viscoplasticity  to  be  reported  in  the  open  literature.  In  this  previous 
research  it  was  assumed  that  the  thermal  response  of  the  plate  was  independent 
of  mechanical  deformations.  Thus,  it  was  possible  to  obtain  the  heat  transfer 
solution  first  and  supply  the  results  of  this  analysis  to  the  mechanical 
solution.  The  finite  element  method  was  emp'loyed  for  spatial  discretization 
of  both  the  thermal  and  mechanical  analyses,  and  finite  differencing  was  used 
in  time.  Newton  iteration  was  used  on  each  time  step  to  account  for  global 
nonlinearity.  This  nonlinearity  was  induced  in  the  the-mal  analysis  via 
temperature  dependence  of  the  thermal  conductivity  and  radiation  boundary 
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conditions.  In  the  mechanical  analysis  nonlinearity  was  caused  by 
viscoplastic  constitutive  behavior,  which  was  predicted  with  a  constitutive 
model  introduced  by  Bodner  and  Partom  111). 

In  the  current  paper  the  model  is  extended  to  include  geometric 
nonlinearity  caused  by  large  rotations.  Furthermore,  a  viscoplasticity  model 
proposed  by  Walker  (12j  is  incorporated  for  comparison  to  results  predicted  by 
Bodner's  model.  The  nonlinear  plate  formulation,  which  was  not  included  in 
the  authors'  previous  paper,  is  discussed  in  some  detail.  Finally, 
sensitivity  studies  are  conducted  using  the  algorithm  to  determine  the  effects 
of  heating  rate  and  plate  thickness  on  predicted  temperature,  displacements, 
and  stresses. 

MODEL  DEVELOPMENT 

The  procedure  utilized  to  solve  the  current  problem  is  similar  to  that 
used  in  other  simpler  mechanical  field  problems.  Therefore,  the  model  is 
presented  only  in  abbreviated  form  in  this  section.  Further  details  of  this 
development  can  be  found  in  references  1I0|  and  113]. 

A  pivotal  assumption  in  the  analysis  is  that  the  heat  transfer  does  not 
depend  on  the  mechanical  deformations.  It  is  well-known  that  in  dynamic 
problems  this  coupling  may  be  significant  even  when  the  deformations  are 
elastic  [14].  In  the  case  of  inelastic  response  such  as  is  considered  herein, 
this  coupling  may  be  even  more  substantial  [15].  However,  the  inclusion  of 
this  term  in  the  heat  transfer  analysis  adds  considerable  complexity  to  the 
computational  scheme  [15].  Therefore,  since''in  the  examples  considered  herein 
the  external  source  is  several  orders  of  magnitude  larger  than  the  internal 
heat  generation,  the  authors  have  neglected  the  mechanical  coupling  in  the 
heat  transfer  solution.  This  issue  is  planned  to  be  the  thrust  of  future 
research  by  the  authors  on  this  subject. 
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As  shown  in  Fig.  1,  due  to  the  uncoupling  assumption  the  solution 
algorithm  may  be  constructed  in  two  stages  for  each  time  step.  First,  the 
temperature  field  is  evaluated  by  using  finite  elements  in  spatial  coordinates 
and  the  Crank-N i chol son  scheme  in  time.  Nonlinearity  due  to  temperature 
dependent  thermal  conductivity  and  radiation  boundary  conditions  is  accounted 
for  via  Newton  iteration.  The  resulting  temperature  distribution  is  then 
passed  to  the  mechanical  problem  as  input,  and  the  mechanical  solution  is 
obtained  by  using  plate  elements  in  conjunction  with  the  Newmark-Beta  method 
in  time.  Nonlinearity  caused  by  inelastic  deformations  and  large  rotations  is 
accounted  for  by  Newton  iteration.  This  procedure  is  performed  recursively  to 
march  forward  for  the  desired  time  span.  The  solution  scheme  is  described  in 
further  detail  below. 


Thermal  Analysis 

An  axisymmetric  finite  element  model,  developed  to  include  nonlinear 
radiation  boundary  conditions,  is  used  to  construct  the  temperature  field  as  a 
function  of  r  and  z  for  each  time  step.  A  typical  two-dimensional 
axisymmetric  mesh  for  the  thermal  analysis  is  shown  in  Fig.  2. 

The  governing  heat  transfer  equations  are  as  follows: 


li  - ’■  1?^  -  ^  S)  =  ° 


1  n  a 


I  oT  ,  sT  —  ^  /xA  tA\ 

kr  —  n  +kr  —  n  =  qar  +  eor(T^  -  T  )  on  r 

01  I  oZ  Z  i 


(1) 

(2) 


where  T  is  the  temperature,  p  is  the  mass  density,  Cp  is  the  specific  heat 
capacity,  k  is  the  thermal  conductivity,  and  r  and  z  are  cylindrical 
coordinates,  as  shown  in  Fig.  2.  Also,  q  is  the  heat  flux  input,  a  is  the 
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thermal  absorptivity,  e  is  the  thermal  emissivity,  o  is  Boltzman's  constant, 
is  the  reference  temperature  at  which  radiation  is  zero,  and  n^  and  n^  are 
components  of  a  unit  outer  normal  vector,  finally,  o  is  the  interior  of  the 
domain  and  r  is  the  boundary  of  the  domain. 

Equations  (1)  and  (2)  may  be  cast  into  a  Galerkin  finiLe  element 
formulation  |16j.  Since  this  part  of  the  model  exists  in  the  open  literature 
|16j,  it  is  not  covered  in  detail  here. 

The  resulting  element  equations  are  of  the  form 


where 


c.  .T  .  +  a.  -T  .  =  q. 


"ij  ‘J 


3N,  bN,  aN,  aN, 


.  .  =  I  ( — I  rk  — ^  ^  rk  dA  +  J  N .  (N,  T,  )  .rdr 

1  j  ar  az  az  '  e  v  k  k^  j 


(3) 

(Ab) 


q  .  -  '  N  .qard"  +  ^  M  i  d  1' 

^1-1^  -  1  r  e 

‘e  'e 


:^c) 


Also,  are  the  finite  element  shape  functions,  A^  is  the  area  of  the 
element,  and  is  the  boundary  of  the  element.  A  linear  triangular  (three 
node)  axisymmetric  element  is  used  in  the  analysis.  Equation  (3)  for  each 
element  is  assembled  to  give  the  following  global  system  of  ordinary 
differential  equations: 

v' 

[Cl  {T}  +  [Aj  {T}  =  {Q}  (5) 

The  Crank-Nichol son  scheme  is  then  applied  temporally  to  obtain  the 
temperature  field  with  time  1 17). 
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Structural  Analysis 


In  the  structural  analysis.  Von  Kartnan 
plate  bending  motion  118),  and  material 
constitution  are  included  in  the  model, 
introduced  via  the  inelastic  strain  tensor, 
next  section. 

The  strain  components,  c^j,  are  defined 


theory  is  assumed  for  the  thin 
nonlinearity  and  viscoplastic 
The  material  nonlinearity  is 
clj,  which  is  described  in  the 


by 


ij 


0 
1  j 


-t-  2. 


1  J 


-5(0.  . 

2  ^  i,J 
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1  .J 


1.2  (6) 


where  c°j  are  the  midsurface  strain  components,  u^  are 
components,  and  are  the  midsurface  rotation  components, 
relations  are  given  by 


the  displacement 
The  constitutive 


=  D.  ..  (cV. 


+  Zk-.j 


1  T  , 

■  'u  ■  'ij> 


(7) 


where  is  the  stress  tensor,  is  the  elastic  modulus  tensor, 
and  elj  are  the  components  of  the  thermal  strain  tensor.  Utilizing  the  above 
in  a  standard  laminate  scheme  will  result  in 


where 


(A.  .,B.  .,D.  .) 
'  ij  ij’ 
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and  h  is  the  plate  thickness.  Note  that  the  coupling  matrix  [B]  does  not 
disappear  due  to  the  through-thickness  variation  of  elastic  modulus,  D^j, 
which  is  temperature  dependent. 

The  governing  equations  for  the  plate  motion  are  thus  derived  by 
satisfying  the  conservation  of  linear  and  angular  momentum  119]; 
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Integrating  equations  (10)  through  (12)  against  variations  in  the  components 
of  the  displacement  field  will  result  in  the  following  variational  principle. 
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where  N  ,  N  .  and  N  are  components  of  the  plate  midsurface  forces  per  unit 

A  y  xy 

length,  u  and  v  are  the  in-plane  displacement  components,  w  is  the  out-of¬ 
plane  displacement  component,  and  p  is  the  out-of-plane  traction. 

Incrementing  the  field  variables,  neglecting  the  third  and  higher  order 
terms  of  the  displacement  increment,  and  applying  finite  element 
di scretization  results  in 


IM]  +  |K1  (au)  =  -  (F^)  -  iF^l 


where  [M]  is  the  mass  matrix  and 


[K1  =  (K^J 


A1  so. 


K,]  =  I  J  IB.,  )  dV 

*-  n=l  V  L  e 


1^1  I  =  I  i  ““e 


n=l  V. 


where  (B[_l  and  IBf^[_l  are  the  linear  and  nonlinear  strain  deformation  mapping 
matrices,  respectively.  Also,  [R}  is  the  external  load  vector  and 
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(18) 


0 

=  I  J  fB|  + 

^  n=l  V  '- 
e 

1aD|  {t°  +  2k  -  -  e^})  dVg 

n 

(^}  =  I  I  IB.,  +  Ae^  + 

n=l  V  ^‘- 
e 

(ADI  (e°  +  Zk  -  -  c^)  dVg  (19) 

and  Hg  is  the  number  of  nodes  per  element. 

Solving  equation  (14)  by  the  (iewmark  integration  scheme  will  give  the 
first  approximation  of  {AUj  at  time  t+At.  The  Newton  iteration  method  will 
give  convergence  to  the  nonlinear  solution  (10,20].  The  authors  are  currently 
utilizing  a  three-node  plate  element  with  five  degrees  of  freedom  per  node,  as 
described  in  references  |13|  and  |21).  We  have  not  employe  an  axisymmetric 
plate  element  because  we  prefer  to  maintain  the  flexibility  to  solve 
nonaxisyrmietric  problems. 

Thermomechanical  Constitutive  Models 

In  order  to  prescribe  the  forcing  functions  [Fj]  and  [F2}  defined  in 
equations  (18)  and  (19)  it  is  necessary  to  determine  the  inelastic  strain 
increment,  This  is  accomplished  by  integration  of  the  selected 

viscoplastic  constitutive  model.  The  authors  are  currently  using  Walker's 
model  112],  as  well  as  the  anisotropic  hardening  form  of  Bodner's  model 
[22].  These  models  are  compared  critically  in  reference  [23]. 

Bodner's  model  assumes 


■I 


(20) 
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where 


is  the  deviatoric  stress  tensor,  whereas  Walker's  model  proposes 


•  I 


(21) 


where 

).  =  .  (o..,  aj,  a^..)  (22) 

and  is  the  drag  stress,  and  is  the  back  stress  tensor.  The  above  is 
supplemented  by  an  additional  set  of  evolution  laws  of  the  form 


£)..  (o,  a,  T) 
ij  '  k£,*  k£’  ' 


(23) 


Equations  (20)  and  (23)  are  typically  numerically  stiff,  so  that 
numerical  integration  to  obtain  is  not  straightforward  124).  Bodner's 
model  is  currently  being  integrated  using  Euler's  forward  method,  whereas 
Walker's  model  is  integrated  using  Euler's  backward  method  |25|.  Both  models 
are  subincremented  within  each  integration  point  on  each  time  step  in  order  to 
produce  accurate  values  for  on  each  global  time  increment. 


Example  Problems 

To  demonstrate  the  use  of  the  algorithm,  an  isotropic  circular  plate  with 
thickness  h=0.042  in  and  radius  r=10  in  is  selected.  The  material  used  is 
B1900+Hf  which  is  a  nickel-based  superal'loy  commonly  used  in  hot  gas 
turbines.  Material  constants  for  Bodner's  and  Walker's  models  are  shown  in 
Tables  1  and  2,  respectively.  The  thermal  boundary  conditions  are  of 
radiation  type  with  reference  temperature  T[^=0°F  on  all  boundaries  including 
the  plate  edge  at  r=10  in.  An  instantaneous  heat  flux  is  assumed  to  be  evenly 
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distributed  over  a  spot  of  radius  0.5  in  at  the  center  of  the  plate.  All 
components  of  the  displacements  and  rotations  are  assumed  to  be  zero  at  the 
plate  edge  r=10  in. 

The  finite  element  mesh  diagrams  for  the  thermal  and  structural  analyses 
are  shown  in  Fig.  2.  Although  a  fairly  coarse  mesh  was  utilized,  more  refined 
meshes  showed  that  the  global  thermal  and  structural  responses  were  not 
significantly  compromised  by  using  the  current  mesh.  The  number  of  global 
degrees  of  freedom  was  held  to  a  minimum  because  the  nonlinear  examples 
presented  herein  required  approximately  15  hours  of  CPU  time  on  a  VAX/8800. 
This  is  due  to  the  fact  that  the  global  nonlinearity  necessitated  the  use  of 
very  small  time  steps.  In  the  examples  shown  herein,  the  maximum  allowable 
time  step  was  0.00001  sec,  thus  resulting  in  a  total  of  1500  time  steps. 

Figure  3  shows  the  predicted  transverse  deflection  at  the  center  of  the 

■p 

plate  for  a  heat  flux  of  60  Btu/in  -sec.  Shown  in  the  figure  are  both  the 
geometrically  linear  and  nonlinear  results  obtained  by  assuming  the  plate  is 
linear  elastic.  The  transverse  deflection  is  significantly  decreased  by  the 
inclusion  of  the  nonlinear  terms  in  equation  (6),  which  demonstrates  the 
importance  of  including  geometric  nonlinearity  in  this  model.  Figure  4  shows 
the  through-thickness  temperature  variation  at  the  center  of  the  plate  for 
various  times.  It  is  apparent  that  a  rather  large  temperature  gradient  occurs 
at  short  times,  and  this  through-thickness  variation  is  the  primary  source  of 
bending  in  the  structure.  Figure  5  depicts  the  temperature  as  a  function  of 
time  on  the  plate  upper  surface  both  at  the 'center  of  the  plate  and  at  r=0.5 
in.  From  this  figure  it  is  apparent  that  a  large  in-plate  temperature 
gradientoccurs  near  the  center  of  the  plate  at  short  times.  This  gradient  is 
the  chief  contributor  to  the  extension  of  the  structure. 
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Figures  6  through  10  are  the  results  for  an  instantaneous  heat  flux  of  90 
2 

Btu/in  -sec  including  both  geometric  nonlinearity  and  the  viscoplasticity 
models.  As  shown  in  Fig.  6,  when  the  inelastic  material  response  is  included, 
the  deflection  of  tiie  center  of  the  plate  is  reduced  by  only  a  small  amount. 
However,  as  Fig.  7  shows,  the  radial  stress  history  near  the  center  of  the 
plate  (r=0.11  in)  and  on  the  top  surface  is  significantly  reduced  by  the 
accumulated  inelastic  strain.  Thus,  material  inelasticity  significantly 
affects  the  predicted  stresses  in  the  analysis.  Figure  8  shows  the  stress 
history  predicted  by  Walker's  model  at  various  positions  on  the  top  surface  of 
the  plate.  Note  that  within  the  radius  of  the  input  heat  flux  the  radial 
stress  attains  a  maximum  near  the  ultimate  uniaxial  strength  of  the  material 
and  then  decreases,  due  to  the  large  thermal  gradient  produced  through  the 
thicxness.  Figures  9  and  10  show  the  radial  stress  and  hoop  stress 
distributions  at  time  t=0.01  sec  on  the  top  surface  of  the  plate.  Both 
components  of  stress  are  significantly  reduced  near  the  center  of  the  plate  by 
the  accumulated  inelastic  strain.  Figure  11  shows  the  in-plane  deformation 
and  Fig.  12  shows  the  accumulated  inelastic  strain  at  time  t=0.01  sec.  The 
inelastic  strain  is  essentially  zero  for  r>l  in,  implying  that  outside  this 
range  the  plate  response  is  elastic.  Also,  Walker's  model  tends  to  accumulate 
more  inelastic  response  than  Bodner  and  Partom's,  thus  predicting  slightly 
lower  stress. 

Figures  13  through  16  show  the  results  for  plates  with  varying  thickness 
and  identical  external  heating.  Figure  13  shows  the  deflection  history  at  the 
center  of  the  plate  as  a  function  of  plate  thickness.  The  center  deflection 
rate  increases  significantly  with  decreasing  plate  thickness.  Figure  14  shows 
the  radial  stress  at  r=0.11  in  on  the  top  surface  of  the  plate  for  various 
plate  thicknesses.  Although  the  radial  stress  peaks  earlier  with  increasing 
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plate  thickness  the  value  of  this  maximum  does  not  change  dramatically. 
Figure  15  shows  the  accumulated  inelastic  strain  at  time  t=0.01  sec  for 
various  plate  thicknesses.  Figure  16  shows  the  temperature  rise  as  a  function 
of  time  at  the  center  of  the  plate  on  the  upper  surface.  It  can  be  seen  that 
extremely  large  heat  fluxes  such  as  those  considered  herein  will  eventually 
cause  localized  ablation  wfiich  cannot  be  handled  by  the  current  model. 

Figures  17  through  20  show  the  results  for  plates  with  fixed  thickness 
and  various  external  heating  rates.  Figure  17  shows  that  increasing  the  heat 
flux  produces  modest  changes  in  the  center  deflection  rate.  Figure  18  shows 
the  radial  stress  history  at  r=0.11  in.  Although  slightly  larger  compressive 
peak  stress  is  predicted,  the  time  at  which  this  peak  occurs  increases 
substantially  with  decreased  heating  rate.  Fig.  19  shows  the  radial  stress 
distribution  at  time  t=0.01  sec  for  different  external  heating  rates. 
Finally,  Fig.  20  depicts  the  temperature  rise  on  the  upper  surface  at  tiie 
center  of  the  plate  for  various  heating  rates.  As  expected,  the  temperature 
rise  increases  substantially  with  heating  rale. 


s 
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CONCLUSION 


In  this  paper,  the  author's  previous  research  [lO]  has  been  extended  to 
include  large  deflection  plate  theory  and  Walker's  viscoplasticity  model. 
Studies  have  been  carried  out  to  investigate  the  structural  response  of  a  thin 
plate  subjected  to  rapid  eAiernal  heating.  Several  conclusions  can  be  made 
from  the  current  results.  First,  due  to  the  thermally  induced  large  stresses, 
the  transverse  deflection  of  the  thin  plate  is  significantly  decreased  if 
large  rotations  are  considered.  This  indicates  the  importance  of  including 
nonlinear  geometric  effects  in  the  current  model.  Second,  the  comparison 
studies  show  that  while  the  two  viscoplastic  constitutive  models  predict 
approximately  equivalent  effects  on  the  structural  response  the  affects  of 
viscoplasticity  are  substantial  and  therefore  cannot  be  neglected.  However, 
Walker's  model  tends  to  accumulate  larger  inelastic  strain,  which  in  turn 
predicts  lower  stress  in  the  structural  response. 

Results  indicate  that  the  deflection  rates  are  increased  with  decreasing 
plate  thickness.  Furthermore,  the  model  accumulates  greater  inelastic  strains 
with  decreasing  thickness.  Finally,  when  the  plate  thickness  is  not  changed, 
the  displacement  increases  with  increasing  heating  rate.  In  this  case,  only 
slightly  greater  peak  compressive  stresses  are  predicted,  but  the  peak  stress 
occurs  more  rapidly  with  increasing  heating  rate. 
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Table  1.  Bodner's  Constants  for  B1900+Hf  126] 


Temperature- Independent  Constants 
=  .270  Mpa'^ 

M2,  =  1.52  Mpa”^ 

=  0.0 

=  3000  Mpa 
Zj  =  1150  Mpa 
ri  =  r2  =  2 
Dq  =  1  X  lO^sec”^ 

Temperature-Dependent  Constants 


Temp  (C) 

n 

7601.055 

2700 

0 

2700 

871  1.03 

2400 

.0055 

2400 

9820.850 

1900 

.02 

1900 

10930.70 

1200 

.25 

1200 

Z2=Z^(Mpa) 
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Table  2.  Walker's  Constants  for  B1900+Hf  [26| 


Temp  (C) 

D{Mpa)) 

V 

Kj(Mpa) 

6(Mpa) 

ri2(Mpa) 

^'3 

21 

1.900E5 

.322 

12.4 

1.73E11 

2.41E6 

4794 

427 

1.900E5 

.328 

12.4 

1.73E11 

2.41E6 

4794 

538 

1.900E5 

.331 

12.4 

1.73E11 

2.41E6 

4794 

649 

1.800E5 

.334 

12.4 

3.862E10 

8.27E5 

1714 

760 

1.655E5 

.339 

13.8 

2.55E10 

8.27E5 

1880 

871 

1.438E5 

.324 

16.6 

5.50E11 

2.36E5 

621.1 

982 

1.249E5 

.351 

13.8 

4.20E10 

9.65E4 

400 

iC93 

1.161E5 

.351 

9. 

5.57E9 

2.36E4 

278.7 

Temp  (C) 

''4 

'■5 

% 

"9 

'^10 

nj 1 (Mpa) 

21 

0 

.3117 

0 

11.87 

0 

4.7E3 

427 

0 

.3117 

0 

11.87 

0 

4.7E3 

538 

0 

.3117 

0 

11.87 

0 

4.7E3 

649 

0 

.3117 

0 

16.64 

0 

4.7E3 

760 

0 

.311  / 

c 

19.83 

2.44E-3 

4.7E3 

871 

0 

.3117 

8.73E-4 

59.33 

2.44E-3 

9.65E2 

982 

0 

.3117 

4.29E-4 

136. 

2.44E-3 

0 

1093 

0 

.3117 

4.83E-2 

136. 

2.44E-3 

0 
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Fig. 5  Temperoture  History  at  Center  of  Upper  Surface 
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Fig. 14  Radial  Stress  History  at  r=0.11  in. 
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A  co«>p**enens  ive  experipe«i4l  p»'oortiri  has 
oeen  oevelopeO  to  Slu0>  tne  transient  response  of 
a  viscoelastic  plate  sucrecteo  to  rapio  heat 
input.  Tne  expe»-imenis  consisted  of  irradiating 
a  clanpec  plate  of  hestelloy  X  witn  a  high  energy 
IS  Kw  c.*».  laser,  with  flisplace»eni  tf^ 

leavoeT'aiure  serving  as  aeasured  data.  a 
tneo^encel  »oOfl  is  also  oeing  developed  fo'r 
co«oarison  to  the  experimental  results.  Tne 
vooel  assumes  one-way  inermomecnanical  coupling; 
mat  IS,  it  is  assumed  that  tne  temperature  field 
is  inoepenoent  of  oeformation  out  not  vice 
versa.  Tne  tnenxal  analysis  is  nonlinear  via 
radiation  Dounoar^  Conditions  and  temperature 
oeoenaen:  thermal  conoociivities.  The  structural 
analysis  includes  geometric  nonlinearity  and 
material  viscoplasticily.  Ine  experimental 
proceoure  and  an  aooreviaied  development  of  tne 
model  are  descriped  in  mis  paper. 

Nomenclature 


-  area  of  two-dimensional  finite  element 

It,  1  -  matrix  relating  linear  components  of 

Strain  and  oelonoation 

!tj^.  1  -  matrix  relating  nonlinear  components 

of  strain  and  deformation 

Cp  -  coefficient  of  specific  heat  at 

constant  pressure 

-  elastic  modulus  matrix 

f^  -  components  Of  oody  force  per  unit  mass 

k  •  thermal  conductivity 

•  comoonents  Of  finite  element  mess 
matrix 

r.r*  •  components  of  unit  outer  normal  vector 
in  Cylindrical  coordinates 

N,  K  -  components  of  plate  miosurfacc  forces 

and. moments  per  unit  length 

•'i  -  finite  element  shape  functions 

d  -  heal  flux  Input 
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-  radial  coordinate  direction 

-  finite  element  external  force 
components 

•  time 

-  temperature 

-  reference  temperature  at  which 
radiation  is  yerp 

-  components  of  plate  displacement 

-  components  of  displacement  in 
Cartesian  coordinates 

-  axial  Coordinate  direction 

-  thermal  aPsorptivity 

-  drag  stress 

-  back  stress  tensor 

-  components  Of  plate  rotation 

-  PouAoary  of  a  domain 

•  boundary  of  a  finite  element 

•  thermal  emissivlty 

•  Strain  tensor 

-  plate  mid-surface  rotations 

-  elastic  strain  tensor 

-  Inelastic  strain  tensor 

•  tnermal  strain  ensor 


-  nodal  componc.-ts  of  plate  rotation 

*  j  / 


X  -  normality  coefficient  for  inelastic 

strain  rate.. 

e  -  mass  density 

c  -  boltyman's  constant 

-  stress  tensor 
V 
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•  •cvUXoric  ilrcu  tensor 


•ij 

6  •  interior  of  I 

Introouct  Ion 


The  transient  response  en  Aerospace 
scrvCti/re  eaposetf  tc  rspic  ntst^hf  be  vast)/ 
different  fro*  that  produced  under  stead/  stale 
conditions,  for  applications  such  as  9as-turbine 
cooDustors.  rockei  n027)es,  space  structures,  ar^ 
me  like.  larae  lenperature  Gradients  arc 
coavaon)/  induced  m  a  ver/  snort  period  of 
tiae.  However,  due  to  tne  conp)ex  and  severe 
envirorment.  )Ut)e  eaperiaenia)  data  are 
ava1)aP)c  for  use  in  developing  nev  inetnods  of 
ani»/si5  Where  aateria)  inetasticU/  occurs  in 
me  structure.  Inerefore.  an  experiaenta) 
proora*  was  esiablishtd  to  invest loate  tne 
response  of  e  avtanic  plate  undergoing  rapid 
heating  with  tetsperau/re  ane  disp)ace«eni  as  the 
primary  var>ao)e5  of  imeresi.  Ine  parameters 
varied  in  the  ekperimehi  included:  1)  p)aie 
tMckness;  2)  heat  f)ux;  and  3)  duration  of  heat 
fiux. 


Thermal  Mairsli 

^  aklSMMetric  finite  elMent  model, 
developed  to  >nc1iide  nonlinear  radiation  boundary 
conditiom.  is  used  to  construct  me  Icmoeraiurc 
field  as  a  function  of  r  and  2  for  each  time 
step.  A  typical  t««o>dimens1oni  t  axis/mmetric 
heat  transfer  meSh  is  Shown  in  fig.  2. 

The  governing  heat  transfer  eouatlons  are  as 
foDows: 

'V  ^  - -h-'"  ‘  '!> 

kr  ^  '‘j  ■  0®^  ^  -  ’T*)  Oh  T  (?) 

Couattons  (])  and  (?)  ma/  be  cast  into  a 
Ga)erain  finite  e)eaent  formulation^.  Since  tMs 
pan  of  me  moot)  exists  in  me  open 
iiteraiurei.  It  is  not  covered  m  detail  nere. 


Tne  resulting  elemen:  eouations  are  of  the 

fonr. 
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A  solution  aloorltha;  has  Peen  developed  b/ 
the  tumors  for  predicting  tne  transient  response 
of  plates  When  sub.ieciec  to  rapid  neat  input. 
wnDe  me  mode)  itse)f  is  not  me  subject  of  mis 
paper.  It  IS  presented  in  aborevieled  form  and  is 
used  to  compart  theory  lo  exoerimen:.  Tne 
solution  algorithm  maxes  use  of  ontowa/ 
rnermomechanica)  coupling  (n  is  assuMd  mat  me 
temperature  fielc  u  independent  of  oef ormatlon) . 
meraoviscoplasiic  responst,  and  oeomeiric 
nonl inearlties.  Tne  finite  element  method  is 
•'tec  v  p  tw?  tttr.  erocesi  to  compute  me  plete 
response.  Tne  thermal  analysis  is  performed 
first,  taking  into  account  such  things  as 
nonlinear  effects  oue  radiation  neat  exchange  and 
lemperaiort  dependent  material  properties.  Tne 
results  are  men  used  as  input  to  the  structural 
anal/sis  wnere  nonlinear  material  response,  using 
faoflner's»  and  Walker’s’  constitutive  models  and 
large  roieiional  plate  theory  are  used  to  compute 
plate  deformation  as  well  as  stress  and  strain. 
Both  me  experimental  and  modelling  programs  are 
describee  in  tne  following  sections'. 

fiQdel  Develpornent 


Tne  model  is  p-esented  in  abbreviated  fonr, 
in  this  section,  former  details  can  be  found  in 
references  5  end  4.  As  snown  in  tne  flowchart  in 
fig.  I,  the  solution  algorithm  is  constructed  in 
two  stages:  tne  meroal  analysis  and  me 
structural  analysis.  On  a  given  time  step,  me 
tnermal  loads  ere  evaluated.  Tnsn  the 
temperature  field  is  solved  0/  tne  finite  element 
metnod.  Using  this  temperature  field,  along  with 
tne  inelastic  strain  increment  evaluated  from  me 
previous  time  step,  yields  me  thermal  strain. 
Tne  themal  strain  results  In  an  unbalanced  load 
from  which  the  oefonmation  field  Is 
eporoximeted.  An  iterative  procedure  is  utHijed 
to  bring  the  solution  to  convergence  for  a  oiven 
time  step.  Tne  solution  Is  described  in  former 
detail  below. 
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A  linear  triangular  axisyim»etr1c  element  is  used 
In  me  anal/sn.  louetior  (3)  for  each  element 
Is  assembled  to  give  tne  following  glopal  system 
of  ordinary  cifferenilal  eouatlons: 


(cl  <T)  .  iAj  (7J  -  iO)  (5) 

Tne  Crans.Hicholson  seneme  Is  then  applied 
temoorally  lo  obtain  mt  temperature  field  with 
time*. 


StruCturt*  Anfelytis 


In  me  strutloral  analysis,  Von  Aarman 
theory  is  assumed  for  the  tnin  plate  bending 
motion',  anc  materi.  nonlinearity  and 
viscoplasTic  constiiutio;  are  incluoec  in  the 
mooel.  Int  material  nur. 'neeriiy  i;  introduced 
vie  the  inelastic  strain  tensor,  it  whicn  is 
described  in  me  next  section,  ^ 

The  strain  comoonentJ  are  defined  as 

Sj  ■  'ij  *  '“■ij  *  I  ’  "j.i’ 


- 1  ’-j  *  1-2  (6) 


Utllliinc  tne  above  In  a  standard  laminate  scheme 
will  result 
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If 
I  B 


)  C2  (8) 

’? 

Nove  lh«l  coupling  »*lrtx  III  ooes  noX 

due  to  the  tnroug^  thickness  vArietlon 
0^  elestic  mooulus,  unich  is  tenper«ture 

oepenoent.  ^  ^ 

Ine  opverninc  eouAtions  for  the  pleie  lotion 
«**e  inus  oerivec  P>  setufyinc  ine  conser^ition 
of  linear  tnc  anouler  noweniuff.  j £} : 

N  «  K  •  OU  *  D  ^01 

*.*  /*•>  .11 


'  ''j./  ‘  ‘“'.ti  •  “  no) 

-  ‘j  - 

^2^  V  •aw  *0  ni\ 

»  ,*k  *y  ,*y  y  ,yy  **  ♦It 

Inteoranno  eauaiions  (9)  tnrooph  (11)  epeinst 
ve»'i«tions  in  ine  cofliponents  of  me  fluoleceineni 
fielc  will  result  in  tne  follcK-'ing  verietionel 
princ'sle. 


tu  1 1'  •  I.  • 

A  ,t 


•  at-  ,)  0*0y» 


j  *»  K  •  OxOy* 

A  *y.*  y.y  .tt'  ^ 

J  tw  (iLjjVj  .  Ej.(c®  -cj^. 
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iV>  'Jj  J  lit) 

Also,  (R)  Is  tnc  eKterittl  lo*4  vector  entf 


Ifji  •  I*  J  IE.l^(t.)  -  lDl'**^*t*  ♦  At*^) 
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Sowing  eouAtion  (13)  by  the  Neumark 
inteoratioh  senene  will  give  the  first 
APprOfinatloh  of  uul  at  ti«e  wai.  Ine  heMion 
kapnsoft  UerAtion  neihoc  uill  give  convergence  to 
the  nonlinear  solution*,'. 

Iheriiomechanical  Constitutive  HogeU 

Jn  oroer  to  prescribe  the  forcing  functions 
(r^)  *ne  (r^)  oeflnec  m  eouatlDns  (17)  ano  (16) 
It  is  necessary  to  oetermine  me  inelastic  strain 
increwent,  atf,.  1ms  is  accomplisneo  by 
integration  ^tf  the  selected  viscoplastic 
constitutive  ittoel.  We  authors  are  currently 
using  uaUer's  aooel',  as  well  as  the 
anisotropic  haroeninc  for»  of  bodner's 
aiooel*^.  Inese  BOOeU  are  coopared  critically  in 
reference  11. 

booner’s  nodel  assumes 

(IS) 

wnereas  welaer's  ®ooel  proposes 

‘Ij  *  ‘‘‘ij  - 


^  ^  ^3* 

end  me  above  is  supplemented  by  an  additional 
set  of  evolution  laws  of  me  foriD 


""'j  ■  <‘icf  ■“«> 


''»'“,xx  *  ^".xy  *  V.yy  ■  "‘'.tt'  ° 

Incrementing  tne  field  variables,  neglecting 
me  third  and  nigher  order  lenns  of  the 
displacement  increment,  and  epplymc  finite 
element  discretiyation  results  in 

.  UJ(.U)  -  (R)'"**'.  (Fj)  .  (Fj)  (13) 
wtiere  Ik)  1s  the  ness  wtrix  ena 
ICJ  •  |r.,J  .  ii:^]  ,l<) 


"=l'  -  I  J  .V 

n*  1  V  ^  ^  ® 


Icuetions  (19)  and  (22)  are  typically 
numerically  .stiff ,  so  mat  numerical  inteoration 
to  Obtain  ij^  is  not  $ir*ightforward»5.  booner's 
Rodel  ij  currently  being  .nieoreicC  using  Culcr's 
forward  rnemoc,  wnerea  walker’s  model  is 
integrated  using  Euler’s  .•JCkward  methodJ-’.  Both 
models  are  suoincremented  on  each  time  stw  in 
order  to  proouce  accurate  values  for  at^.^on 
each  time  increment.  *  ^ 

EKPfrjmental  Prooram 

In  order  to  verify  analytical  models 
currently  being  constructed,  an  experimental 
program  has  been  developed  to  investigate  the 
transient  response  of  a  viscoplastic  plate 
subjected  to  rapid  heat  input.  Of  particular 
interest  is  the  measurement  of  displacement  and 
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lcv>erjiare  f>e)es  for  t  r^ci^n^'isr  p)4t* 

%PtC\mer)  wnoer9otf«( 

expert^nt«*l  pro9r*»  fs  PW^omS  Into  v«o  pn*sct: 
})  t  )ow  energy  |<lDOw/^n‘^)  tnput  pnist  orin; 
conouctce  leitv  **>’  Un>*ertUy;  «n0  ?)  i  h>9h 
energy  (>lnw/ir‘)  input  piwse  being  performed  «t 
tne  >^tr  force  wnent  Aeroneuttcs  Loooretoriet 
using  the  L»*<£L  1  ( L4ser-h«roenefl  H»ieri4U 

[vA^utlior.  lADOretpry)  f«c1^Uy. 

Lo>  {nero»  Testinc 

The  lo*  energy.  sewi-rtcitf  heating 
expermenis  mane  use  o<  «  oane  o<  ouart?  Unts. 
to  irr*aietf  me  plate  specmens.  Ine  purpose  of 
these  ejiperi«ents  is  twc  folc:  first,  mey  serve 
4S  t  SiiihjUtiori  of  me  moh  energy  experifienl* ; 
and  seconc,  they  proouce  ‘e«oeraiure  and 
cisplacenenl  fielcs  wnicn  Change  at  a  fcuch  slower 
rate.  Therefore,  ve'iuaole  mformatioh  with 
repere  10  e*peciec  e*per mental  results  can  oe 
oeinec,  in  aocmon  to  verifying  inooel 
preCiciions  unoer  less  severe  conOilions. 

Ihe  heat  source  for  these  exoe*“ie»ents  are 
toOular  Quartz  larps  with  tungsten  emitters.  A 
circular  neat  target  zone  is  obtamec  by 
shielcmg  excess  energy  froff>  me  specimen  via  a 
water  coolec  orifice  blate,  wtiicn  is  fuieC  with 
a  variable  size  aperture.  inu  apparatus  is 
attacnec  PircCtly  to  me  specimen  Support  fixture 
as  Shewn  in  Fic.  j.  The  support  fixture  imposes 
ClaKii-1ik.e  bounoary  concuioni  along  each  eoge  of 
a  rectangular  specimen  (details  on  tne  support 
fixture  are  discussed  later  m  me  paper), 
insulation  shields  enclose  pom  sides  of  me 
specimer.,  in  oroer  tc  reduce  the  amount  of  neat 
transfer  oue  10  free  convection. 

Displacement  and  temperature  oaia  are 
o&tamec  usme  a  PC  oasec  date  acouisition 
syster,  «unich  also  serve:  as  me  ouert;  lamp 
contrcHer.  A  single  OC-opere:ec  IVDT,  wnicn  can 
be  oositionec  a:  various  racial  location:,  is 
usee  ;c  metsure  displacements,  A  coupling  tang 
enables  dual  measurement  of  displacement  with  me 
LVD*  and  a  precision  dial  indicator,  «easuremen: 
cf  tne  in-plane  temperature  field  and  transverse 
temperature  orafiients  arc  maoc  using  K-type  30 
cage  mermocouples  intrinsically  mountec  tc  m« 
plate  speCimer.  Tne  specimens  for  this  Phase  of 
testing  were  mace  of  6061-T6.  alutninun..  This 
material  was  selectee  because  1:  exhibits 
inelastic  oeforraation  at  relatively  low 
tenpe^atu-es  (^SOCT). 

Testinc  to  oate  has  not  yieloed  a  hign 
enoucr.  input  energy  tc  produce  a  significant 
specimen  temperature  rise.  Jn  addition,  me 
current  aperture  iysieir-  has  failed  to  generate  a 
well  defined  tares:  zone.  Tnerefore,  me  orifice 
plate  is  being  refitted  with  a  lens  syster  me: 
will  serve  tc  Dcm  concense  available  neat  energy- 
end  focus  it  on  me  target  spc:. 

Hipn  Enerev  Testinc 

The  high  energy,  rapid  heating  expcrimenis 
mace  use  of  an  electric  discharge  15  icw 
continuous-wave  carbon  dioxide  laser  ooeretinc  at 
a  wave  length  of  10.6  ym  with  a  flat-top  beam 
profile,  tc  Irradiate  me  plate  specimens.  The 
test  apparatus  was  placed  in  e  nitrogen  flood  box 
during  irradiation  in  order  to  prevent  oxidation 


pf  the  toocimn  %un  ^ct  while  It  wes  being 

heetec.  In  additioh,  ihstrumentat  ion  tor 
■cesurimg  temperature  arid  disolMcameni  fields,  as 
1^11  as  hsgh  speed  photographic  documental  ion 
i^rt  used  and  will  pe  discussed  ii>  detail  below. 

A  ipecimen  support  fixture  was  designed  to 
imoote  clamp-liae  boundary  conditiohs  along  each 
edge  of  *  rectangular  plate  specimen.  Trie 
superstructure  of  me  fixture  was  fabricated  from 
6061-16  aluminor  to  support  the  specimen  in  the 
vertical  position,  as  tnown  in  Figs.  4a  ane  4t. 
Ah  insert  aiaoe  pf  3C>d  stainless  steel  serves  both 
as  a  swpprri  stiffener  and  water  yacxet.  a  total 
waste  water  cooling  sysie*  is  used  tc  proviot  a 
uhiiore  aftC  constant  plate  oounoery 
lemoeraturt.  Inenaocouples  were  attached  to  earn 
half  of  me  insert  at  me  mio-po: -t  lor., 
approximately  0.5  m  away  fro*  me  specimer.,  to 
recorc  any  temperature  variations. 

The  two  piece  (with  me  inserts  m  place) 
picture  frame  style  fixture  securely  clamps  the 
specimen  m  position,  using  Ti  socxel  neac  cap 
screws  toroued  to  1?D  m-lbs.  An  ofi-set, 
staggered  screw  pattern  is  used  to  ensure  a 
uniTora.  cUmoihc  zone,  in  addition,  the  fixture 
IS  indexed  with  naroened  dowel  pins  so  mat 
•  Ugna»ent  between  the  fixivrt  and  specimen 
remains  consistent  froa.  test  to  test. 

The  material  selected  for  mi$  reseerett  is 
Hastelloy  X.  tiastelloy  2  is  a  nicxel-ChromiuB- 
iron-molypoenua.  alloy  mat  Possesses  a 
comomatior.  of  oxidation  resistance  and  high 
strength  at  temperatures  m  excess  Of  ??0CT.  it 
was  selected  for  this  test  program  pr)»ar11y  for 
us  high  temoerature  characteristics,  but  also 
because  It  is  widely  used  in  the  research 
comiunily  in  conjunct  ion  wUn  high  temperature 
constitutive  eiodelinc  and  testing. 

The  hastelloy  X  material  wat  obtained  in 
plate  fonr  u.  no»t.tnal  inicxnesses  of  J/J6  ir..  }/B 
ir..  and  }/<  in.  and  was  usee  to  fabricate  t  total 
of  16  Specimens  (6  of  each  thickness).  The 
material  was  received  in  an  annealed  confliiioh 
soedfieo  by  ASlK  5536  and  used  without  further 
heat  treatment.  ho  micrograohlc  studies  were 
performer  tc  investigate  the  variations  in  grain 
structure  or  sizt  thtt  existed  between’  me 
different  plate  thicknesses. 

All  specimens  were  machined  to  finished 
dimensions  of  13  x  13  in,  corresponding  to  the 
o-te**  eoge  of  support  fixture.  Therefore,  taking 
mto  account  me  support  fixture  damping  zone, 
the  effective  plate  dimensions  were  reduced  to  10 
X  )0  In.  Each  specimen  was  hand  sanded  with  /3?0 
c’-it  sane  oaper  and  cles'  ocad  blasted  to  produce 
a  low  luster  finish.  Nv  other  soecial  surface 
treatment,  sucr.  as  anoi;..-n5,  blueing,  d'  the 
like,  was  performed  to  enhance  the  ihermophysical 
properties  of  me  maieri^t'!. 

I 

An  Integrated  Ins^rjvnvnTation  package  was 
used  to  simultaneously  measure  me  displacement 
and  tenocralure  fields  of  a  plate  specimen 
undergoing  lescr  irradiation.  The  primary 
instrumentation  .ncluoec:  })  LVDT's  (Linear 
Variable  Differential  Transformers)  for  measuring 
displacement;  2)  thermocouples  for  measuring 
tewoeraturc;  3)  a  radiant  pyrometer  for  measuring 
surface  brightness  temperature;  and  <)  strain 
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frequency.  ^  12  til.  ^)9l^  spe«e  MU  %yi\.m 
(c*”''>ec  the  PO*.)  w«s  wvec  lo  cofwe^t  i»»e  •rv«lo9 
output  0^  tr«ntouc«ri  tnc  thermocouples  tc 

#ri  eoutvtienl  tiMry  tore  At  Ati  APDroj^Mte  retc 
Of  1.2  IKZ.  Ohcc  converiec.  the  oaU  were  stored 
Oh  magnetic  tApe  for  supteoueni  concert  tori  ix 
rnginecrtn^  ohUt  Anc  Any  other  poSl- 
processtn^.  Descrtoee  oe1o»  is  i  here  oeiaHeO 
tiitcusston  on  tne  t«ole*enut ton  of  toe  vartout 
pieces  of  >r>tt«v«eniAt  ton. 

loiel  of  21  Ot-ooertitrO  LVC^'s  wert  useC 
ic  lenie  tne  out ^pf-p lane  otsplacements  resutttng 
frof  int  iaser  oePOSMion.  The  outputs  of  ine 
lvDf‘s  were  scalef  vu  mt  oati  svster,  ic  oetect 
CUP  lacement  s  as  smell  as  C.OOOl  in,  ai  t 
puplunec  *ia*taKAr  freooency  response  of  li  h2. 
Tne  LVDVs  were  arranpefl  tr  1  svmmetrtc  pattern 
arounc  me  heal  2one  anC  were  rtcidl)  mountec  to 
a  support  sy’SieV'  wnten  wet  posit'onec  Orectly 
Peh'jnc  me  specimen,- es  snowr  sn  fic.'i  &e  anc 
bt.  Ir  AOdUion,  7  IVOR’S  wfrt  usee  to  montlO’' 
reUtive  movement  oeiween  tne  IVCT  support  sysie<r 
anc  date  f uiurc. 

Heasuremenl  of  me  m-pUne  temperature 
flelc,  mrouQf.-tniceness  leaoeralure  oraoients, 
anC  non-contact  p)aie  surface  temperatures,  were 
meOe  usinc  22  f-lype  30  oape  ineraocouples.  Tne 
tnermocouples  were  concentreicC  malm 
Ctameter  circle  arounc  me  neat  7one  ano  were 
a'-ranpec  m  a  symmetric  pattern  for  reasurmc  me 
m>plane  temoeratufe  fielc.  Tne  tnroupn- 
tMcKness  tefroertture  oraflienis  were  •*easure£ 
usme  thermocouples  posUionefi  at  the  savi^ 
CDOrcinatc  location*,  out  mountec  on  the  front 
anc  oaCi  of  tne  soecimen.  AU  tnerirocouolcs  were 
Intrinsical )y  mountec  lo  me  specimen  via  a 
welfimc  operation  with  me  excection  of  i 
inermocoup'ie* .  T*-ese  thermocouples  were  peac 
weloec  functions  wnich  were  mouniec  approximately 
C.C&C  m  from  me  surface  of  me  plate.  All 
inermocouples  sountec  on  me  front  surface  of  tne 
plate  (the  hea:  sloe)  were  Inconel  60G  sneatec  to 
vftnstanc  :ne  extreme  temperatures,  whereas  the 
thermocouples  mountec  on  tne  oeci  sioe  of  me 
soectffien  were  msulatec  usinc  Moh  teftioereture 
class  Praic.  Tne  thermocouples  were  connecteC  to 
the  oata  system  via  a  25C*f  reference  oven,  which 
for  IMS  test  was  left  open  tc  room 
lemoeraturt.  For  thermocouple  input,  me  oata 
system  was  scaled  to  recorc  voUaoe  cnanoes  on 
the  oroer  of  C.C3  mv,  wMcn  corretooncs  tc  a 
iT>easu-ec  teroeraturt  resolutior  c'  aprroi'ifr.ately 

Thus,  tawmc  irtc  account  tne  A’Si  wire 
error  Specification  anc  tne  aoove  resolution,  a 
maximum  temperature  uncertainty  oeiween  S.3tT 
end  can  oe  expected. 

A  oermantum  radiation  pyrometer  was  used  to 
obtain  relative  measurements  of  me  plate  surface 
priontness  temperature.  Tne  pyrometer  is  a  hion 
soeed  trenscucer,  having  e  peak  spectral  response 
at  a  waveler>ath  of  l.Suff-  and  an  effective 
temperature  range  between  90G*F  end  5400’^  within 
a  target  area  of  eoproximatcly  0.0491  ir,^,  Tne 
pyromeie*’  was  aligned  to  record  temperatures 
within  me  laser  irradiated  soot  diameter  in 
conjunction  with  a  thermocouple.  Tne  output  of 
the  pyrometer  was  fee  into  tne  data  system  for 
use  later  in  oeveloping  en  appropriate  transfe- 
functioh  for  tne  slower  responding  tnemocouples. 


A  36C  t  strafn  9A9e  was  used  tc  mcAWnr  the 
OyriAMlC  respOAAc  of  tiK  stKimen  rcVuUlhg  fro* 
the  rAptc  heAtSng.  The,  strAtn  l^Ape  was  iMwhttC 
pArAliel  with  the  edge  of  the  plate  toecimen 
APprotiwAiely  2.5  lh  off  center.  Since  vlprAtion 
frequencies  And  not  strain  mAOhluioes  were  of 
interest,  tne  Mtput  of  the  strain  gaoe  was 
displAyeo  or.  a  vUlCoroer  strip  Chart  recorder. 

A  IoIaT  of  a  lASer/Siructore  InterACtion 
experiAients  n«vc  been  conoucied  mws  far^ 
representing  AporouMlely  one  third  of  those 
plennec  unoer  this  pnast  of  testing.  All  <  tests 
utilUeC  me  SAOc  l/2t  tn  imc*  plate  specimen 
with  neat  flua  anc  exposure  time  sei^vinc  as  input 
test  variables  (as  snp«.v>  ir  laple  1)  anc 
ciSPlaceAifn:  and  temperAiure  fielos  serving  as 
measured  o'  output  ouan: iiies.  keipre  testing, 
the  laser  oeah  was  maracteruec  m  terms  of  me 
AAOum  of  energy  being  oelivered  or,  target^  in 
addition  to  tne  bear  width  ane  oensity.  large* 
energy  was  measurec  by  reflecting  approx  i».aie  iy 
7S  o'  tne  total  available  energy  into  t  tprpeoo 
calorimeter.  The  bear  brofiU  was  oeienrmec  oy 
Ablation  of  plexiglass  samples  as  snewn  u.  Fig. 
£.  by  measormc  me  plexiglass  burn  patterns, 
both  me  laser  target  area  anc  bean  wnifonmily  in 
the  radial  direction  can  pe  found.  For  mis 
experiment,  the  User  contact  spot  was  found  to 
be  ellipsoidal,  navinc  ma.ior  and  minor  axis 
length  eoual  to  C.5235  m  and  C.5?)6  in, 
respectively  finis  U  a  result  of  me  beam 
sirixiftc  me  specimen  a:  a  2D*  mcioence  ancle  m 
orocT  to  prevent  energy  feedback  inrougn  me 
User  pattern). 

Figures  7  thrwjph  10  Show  typical  output  of 
me  menftocooblet  an/?  iVDT*s  at  oiscrtie 
locations  curing  laser  irradiation  (with  anc 
without  lima  averaging).  In  particuU*-,  figs.  7 
anc  5  Show  raw  oata  (in  engineering  units)  wnere 
mere  n  an  aoparent  electrical  nout  prooler..  A 
power  spectrum  oensity  analysis  mdicatei  mat  $D 
M2  and  its  harmonics,  -ere  present  in  the  oata. 
This  resulted  m  peai-to-peax  venations  of 
Approximately  Aa*f  end  £.0074  in  for  temperature 
Anc  oiSblAcemcnt,  respect ively.  future  tests 
will  incorporate  real  time  analog  filtering 
and/o**  post  test  digital  filtering,  in  oroer  to 
resolve  this  problec.  Figures  £  anc  20  present 
tne  seme  oata  as  discussec  above,  except  mat 
tiR^  averaging  has  been  usee. 

Ko  conment  is  made  about  the  validity  of  the 
oata.  except  to  note  tnt  fpllo-inc.  First,  tne 
typical  temperature  profile  snewn  m  rigs.  7  anc 
£,  indicates  tne  thermocouple  die  not  respond 
during  tne  laser  shot.  Tms  should  no:  be 
interpreted  as  tnermornuple  lac.  For 

intrinsically  mounted  tt.  'mocouples.  one  would 
expect  a  response  time  m  the  msec  time  frame. 
It  IS  presumed  tne:  me  user  firing  seooence 
caused  a  voUge  shift  an**  saturated  the  A  tc  0 
system  low-.  Second,  t  ta'get  rone  temperature  of 
8b0-200C*f  was  generally  encountered,  which  was 
far  below  tne  )5OD-?0X*F  temperature  expected. 
Lastly,  Figs.  9  and  20  snow  a  typical  LVDT  oata 
that  indicates  a  drastic  disolacement  reversal 
flurino  irradietion.  That  is,  the  Specimen 
initially  deforms  in  me  positive  direction 
/towards  the  heat  source)  anC  men  reverse 
itself,  going  through  a  neutral  point  to  assume  a 
negative  displacement.  The  authors  ere  ewaitinc 
Dhotographic  confirmation.  via  high  speed 
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cwrti,  te  let  tett  untnorntru  ecvuili/ 

(ixtl. 

The  pri«*ry  Otiectue  o'  tr,M  'int  P««ie  o' 
leilint  .«s  tc  tJliolilo  i  witH'»  eipenoenlit 
proceaurt  kMC*!  •<oo1e  vt'-Uy  inilytlcil  ■oOeU 
Cu'TfnHy  unoe'  COnSlroCt'on.  1c  thU  enC.  tMt 
Defective  ixs  0000  •«,  “'l''  t"'  e»«Pt<on  0'  t« 
toHo«<n{  orotie*  e'eis  (*MCh  ft  currenl  ly 
oeins  torrtcied);  1)  entrciy'nc  ina/or  fuin;  ti>e 
Uver  eowti-iely  *'<00100  tne  M*iortc  0*11.  BOlh 
wU'  noise  *X  »oU*9e  snl'ls;  2)  tnett  yii  "C 
explicit  iiw?ic*tior.  o'  wncr  irt  lesfr  poxer  wes 
or.  *nc  0"  me  H'oei  jete;  3)  me  nuroaer  f  looc 
DO»  01c  not  provioe  *r  me':  environment  *nc  woe 
pnoQOO'ftpnic  ooctmient*  t  ion  ci'flciil:;  d)  tne 
inei-iocoupte  ee'enenct  yunciion  neeoec  to  oe 
pelo»  eoo«  ie«ipe'*tort;  6)  e  consten;  loe'iure 
setting  on  me  vioec  »ni  tilon  soeeC  c*»er*s  aaoe 
pnotoon«n!C  interpretenon  eif'iculi  *s  tne 
soecliien  heetec  up;  6)  tne  pyeomeie'  flett  •*: 
in»oo'Oori»ifly  scelec;  *nc  7)  tne'e  -es  »Doanen: 
wvoment  o'  tne  iWT  support  suni  auring  l*so' 
deposition. 

Coape^ispn  c*  tc  typf**iwen: 

Lm'Uec  nesuVs  nevt  been  obteinec  me 

fOOe"'  tO’’  cowpenuon  to  me  e*penii>eni*'i  rcsuUs 
cwen  m  me  pi'evious  section.  tipen^nul 
constants  fpr  bopner ' s  »ooel  nave  beer  oDtemec 
tor  hasteUoj  X'^.as  sno-r  in  Table  2.  Tne  e»OdeT 
nas  been  useC  to  predict  tne  response  of  tne 
SQuart  plate  of  l/li  m  thickness  wUh  SMtUlly 
end  ipaoorall>  constant  nett  mput  of  12.2 
Etu,Mr.‘  sec  over  t  C.?S  in  radius  soot  at  tne 
center  o'  tne  plate.  freflictec  lemoereiure  at 
tne  f  ront  center  o'  me  plate  versus  tune  is 
snowr  fo'  tne  first  C.Ofc  sec  m  fic.  T..  ^-nere  it 
car.  be  seer  thtt  tne  plate  is  neatec  rapid  ly  tc  a 
iemoe’‘aiure  o'  about  130C*r.  This  neatinc  is 
sufficient  wO  cause  a  cente'’  disolace'^er:  of 
a2P*‘0« iTr'diely  C.C2  m  {aaout  33*  of^  the  plate 
thickness),  as  shown  in  fic.  12. 

Cisdacenent  field  m  turn  proouces  radial  and 

noop  stresses  near  tne  plate  center  wnicr.  rapidly 
exceed  tne  waterial  yield  point  on  me  front 
surface,  as  snown  in  Fios,  13  and  14.  These 

stresses  ere  mouPec.  m’  Urge  measure.  p.v  me 
plate  oenoinc  wnicr  occurs  near  tne  center  of  tne 
pleifc.  as  sn^n  in  Pic.  IS. 

A  co«oarison  between  Fics.  £  end  11  snows 
that  the  predictec  anc  expe’^iflienttlly  ooserved 
ter.pe-atures  e-e  np:  \r  agreement,  ru-r^nennore. 
Pies.  1C  and  12  oemonstrate  tnat  tne  cisDlacen>ent 
predictions  are  else  not  in  accord  with 

eapc’*iff»eniel  oPse*’vetions.  T^e  authors  hasten  to 
pcint  out  that  DCth  me  experimental  and 

theo’*etical  results  ere  first  passes,  and,  as 
such,  np  ettcitpt  nes  oeen  mcoe  tc  ‘messaoe 
either  result  tc  ff'dtch  me  omer.  Given  me 
complex  nature  cf  Poth  me  eoce  .  anc  exoerimen*. 
the  ffiScoroant  results  are  to  be  expected.  Since 
tnere  are  numerous  sources  for  tne  Cisaoreement, 
our  current  empnasis  is  on  uncovering  me  most 
important  sources  m  such  a  way  as  tc  enhance  me 
,*  f  ■  fsf  tne  mooel. 

briefly,  we  suspect  me  following  sources  of 
disaercemen:  oetween  me  model  and  experiment: 


1)  m  Mve  rea*»»  te  wapect  mat  both  me 
LVm'i  end  inerwxcbup  les  art  giving  spurious 
rewHs  Short  t»of  spans.  i*e  currently 
stuoying  this  problem. 

7)  The  oooel  assumes  mat  me  input  source 
coMences  instantaneously,  i^ereas  me  laser 
output  moicaies  an  approximately  linear  rise 
tiw  of  .003  sec.  ife  are  acxJlfying  me  oooel  to 
account  for  mis. 

3)  Initial  if^Jlcatlons  art  mat  me  input  beat 
source  IS  not  spatially  ncpkweneous .  thus 
prooucmc  a  temperature  field  wnich  is  not  a>i- 
sy<nHeirit. 

4)  Currently,  we  can  heasure  only  me  amount 
of  neat  energy  oemc  Delivered  on  target,  future 
tests  will  try  10  Ptasore  tne  amount  of  energy 
ocme  reflected  by  me  specimen. 

5)  Tne  very  rapid  htat  lf*ul  may  cause  tnera.al 

waves  Which  are  not  predicted  by  ov  par^poUc 

neat  equation.  we  are  thus  considering  various 
nyperoolic  forms  ©f  tne  neat  eouation. 

C)  Tne  co4Bputer  rode  currently  reouircs  one 
C»‘II  o*y  or.  t  VJ.X  B88C)  (3  't*l  7’"'  o.-'X)  70 

P'oauce  7m  firs7  Ci.Ot  lec  of  prtaic7to 

'eipof>st.  we  «rt  *tie«otinc  IP  improvt  7ne 
co«oul*;ion*l  efficiency  of  7ne  K>oe1  so  li.*l 
lonoer  lent.  preCic7ioni  c*n  be  oDl.ineC- 


Condlusipn 

It  Is  obvious  tnat  me  experimental  results 
are  not  in  aoreemeni  with  tne  model  at  thi*  pom, 
in  time.  While  me  authors  cannot  say  with 
cerctir.-.v  wnf.  ine  c*uiet  fo'  cisc'ecsncy  *re.  17 
n  me  o’ene'.l  feel  me  nere  7bet  me  e«oe'iiwn7* 
P'oor.ir'neeos  70  pe  refmec  consioerio’iy  befo't  e 
C'liicel  cono.'iiO'.  of  me  Boot)  me  erof'imen:! 
.ill  yitlc  f'uilful  reiulti.  we  ere  ou"en7  ly 
piAicinc  prooress  m  OOth  flireclloni. 

»ctno.leooe»ien7 

Tnil  eesefCP  ues  toonsoreC  Oy  me  fP'Ct 
Office  pf  icien-.lfic  koieerco  unoer  conl'eci  no. 
f«9t?C-86-l;.001t.  3ne  eulno's  woolo  especit'lv 
line  '.0  tnenk  Dr.  Deor?*  P.  Senoecfy;  me  nis 
lupoort  tecnnlciens  el  AFUAL.  Ineif  cominueC 
efforts  brought  me  experimental  phase  of  mn 
research  to  fruition. 
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Table  1.  Laser  Parameters  for  the  Specified  Tests 


■ 


She: 

Jncioen: 

bear  , 

AumDer 

{nergy  (r.W) 

Area  (cff' 

<9540 

1.4 

J.4 

49S4} 

1.4 

4.2 

1.4 

<S6S0 

7.0 

1.4 

Table  2.  booner-Pertoo.  Model  Materiel  Constants  for  hastelloy  X 


Constant 

Temper.iure 

1.0 

kooio  TeAip 

i.C 

1200T 

1.0 

leoo'f 

0.75 

lC"i 

10:4 

ICii 

10{4 

MPa'* 

(Ksr*) 

0.02 

(1.36--<) 

C.l  (6.e52.P) 

0.) 

(6.69; -<) 

HPi‘‘ 

KPe"^ 

(i:sr2) 

1.6 

(:.2<i-2) 

2.t  (1.652-2) 

2.« 

(1. 652-2) 

(Ksr‘) 

C.OOl 

(6.e5;'-6) 

c.oj  (i.ea-s) 

.0) 

(6.652-5) 

HP, 

(Ksn 

1200 

(:.7ct5) 

1550  (2.2CE25) 

•  :o 

(5.72£) 

MPft 

(«i) 

2000 

(2.9-5) 

2000  (2.9-5) 

2JOO 

(2.925) 

MPe 

1200 

(1.7P55) 

3550  (2.2<E25) 

1  000 

(E.724) 

MPt 

sec"^ 

sec-2 

(«i) 

1200 

(1.7P55) 

BOO  (1.1625) 

500 

(7.252£) 

0 

6.52-7 

6.52-7 

c 

6.52-7 

6.52-7 

C.9E 

0.98 

C.9E 

0.98 

C.96 

e.ge 

0.9S 

0.98 

.  MPt 

(«i) 

2.0726 

(2.00227) 

1.6126  (2.22527) 

1.2726 

(:.9B7l7) 

t)  Support  Structure  -  Slot  Vie*. 


P)  Support  Structure  -  Front  Vie>. 


Fic.  i  rUtf  Soecimer  Support  fiitture 


i)  LVDT  Support  Sunc  -  too  View 


b)  LVDT  Support  Stsnd  -  Sioe  View 


Fic.  6  LVOT  Support  Stend  Positioned  behind  the  Specimen 


»)  before  AbUtion  b)  After  Abletion 

Fic.  6  Plexicless  Samples  Used  to  F'.eesure  the  Leser  beam 
Uidth  and  Density 
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